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Abstract 

The data taken with the ANTARES neutrino telescope from 2007 to 2010, a total live time of 863 days, are used to 
measure the oscillation parameters of atmospheric neutrinos. Muon tracks are reconstructed with energies as low as 
20 GeV. Neutrino oscillations will cause a suppression of vertical upgoing muon neutrinos of such energies crossing 
the Earth. The parameters determining the oscillation of atmospheric neutrinos are extracted by fitting the event rate as 
a function of the ratio of the estimated neutrino energy and reconstructed flight path through the Earth. Measurement 
contours of the oscillation parameters in a two-flavour approximation are derived. Assuming maximal mixing, a mass 
difference of Am? 2 = (3.1 ± 0.9) ■ 10 3 eV 2 is obtained, in good agreement with the world average value. 

Keywords: neutrino oscillations, neutrino telescope, ANTARES 



1. Introduction 

A measurement of the quantum mechanical phenomenon of neutrino oscillations provides important information 
on the mass differences of the neutrino mass eigenstates and their mixing angles. The effect was discovered by 
Super-Kamiokande 0]] on neutrinos produced in the Earth atmosphere. In this paper, the first observation of neutrino 
oscillations by a high energy neutrino telescope is presented. As this measurement addresses a higher neutrino energy 
range, the analysis presented here is complementary to previous measurements. In particular, the selection procedure 
for events, on which the oscillation parameters are measured, is different with respect to a similar analysis at lower 
energies. Whereas the separation of charged current events from v e and neutral current events becomes much easier, 
handling the background from misreconstmcted downgoing muons is more challenging (see Section|6]l. 

The main goal of the ANTARES deep sea neutrino telescope [2f] is the observation of high energy neutrinos 
from non-terrestrial sources. The telescope is optimised for the detection of Cherenkov light induced by the passage 
of upgoing muons produced in charged current interactions of neutrinos at TeV energies, for which the muon can 
traverse completely the equipped detector volume. At lower neutrino energies, both the cross section and the muon 
range are smaller resulting in a decrease of the detection efficiency. For the analysis presented here, the lowest energies 
of detectable muons from neutrino interactions are about 20 GeV. For upgoing atmospheric neutrinos which traverse 
the Earth, this threshold is low enough for the observed flux of induced events to be significantly suppressed by 
neutrino oscillations. By studying the observed upgoing muon rate as a function of the ratio of the reconstructed muon 
energy and zenith angle, constraints on the atmospheric neutrino oscillation parameters are derived. 

The paper is organised as follows: in Section 2 the phenomenology of atmospheric neutrino oscillations is intro- 
duced and the method used to extract the oscillation parameters is discussed. The ANTARES telescope is described 
in Section 3. The data set and the Monte Carlo simulations are explained in Section 4. The methods to reconstruct 
the neutrino direction and energy are discussed in Section 5 followed by details of the event selection in Section 6. A 
discussion on systematic uncertainties is given in Section 7 and the final results are presented in Section 8. 
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2. Neutrino Oscillations 



The survival probability of atmospheric in the framework of three flavour mixing is given as 

, ,/1.27Am 2 L\ , , ,/1.27 Am 2 ZA , , ,/1.27Am 2 L\ 
^ -» Vp) = l-^fl^fsin 2 ! _^j-4|i7 Ml | 2 |C/^l 2 sin 2 l I - 4|I/^| 2 |I/^| 2 sin 2 l 

(1) 

where L is the travel path (in km) of the neutrino through the Earth and E v , its energy (in GeV). U a i is the 3x3 PMNS- 
matrix which describes the mixing between flavour eigenstates v e , v M , v T and mass eigenstates vu vs> ^3 and Am?. = 
|m 2 -m 2 | (in eV 2 ) is the absolute difference of the squares of the masses of the corresponding neutrino mass eigenstates. 

The present analysis is restricted to L < 12800 km and E v > 20 GeV, for which the term sin 2 (1.27Am 2 1 L/£' v ) does 
not exceed 3.7 • 10~ 3 when using Am— from [3]. This term can safely be ignored as well as differences between Am?j 
and Am 2 2 and EquationQ]simplifies to 

, , ,/1.27Am 2 L\ 
P( Vfl -» V/1 ) = 1 -4(1 - \U^\ 2 )\U» 3 \ 2 sm 2 [ ^ 2 -j (2) 

Results could in principle be extracted in terms of \U/a\ 2 and Am 2 2 which are the two oscillation parameters in Equa- 
tion|2] To maintain compatibility with earlier results, a mixing angle sin 2 023 = |E^u3l 2 is defined, ignoring the 2.4% 
deviation from 1 of cos 2 0^ = 0.976 fl. This leads to the usual two-flavour description 

,/1.27Am 2 L\ , ,/ 16200 Am 2 cos 0\ 
I E~) = 1 ~ Sin " 2623 ( " (3) 



For upgoing tracks L is in good approximation related to the zenith angle by L — D ■ cos where D is the Earth 
diameter. The transition probability, P, depends now on only two oscillation parameters, Am 2 2 and sin 2 2023, which 
determine the behaviour for the atmospheric neutrino oscillations. 

With AmL = 2.43 ■ 10~ 3 eV 2 and sin 2 2023 = 1 from one expects the first oscillation maximum, i.e. Piy^ — > 
Vfj) = for vertical upgoing neutrinos (cos = 1) of E v =24 GeV. Muons induced by 24 GeV neutrinos can travel up 
to 120 m in sea water. 

The observed number of events in bin i, N,, of a given variable can be compared to the number MC, of expected 
Monte Carlo events in the same channel 



1 - sin 2 26 n sin 2 



16200 Am 2 2 cos ® ik 



(4) 



where fi, is the number of the background atmospheric muon events in channel ; and the sum gives the number of 
atmospheric neutrino events in channel i weighted by the event dependent oscillation probability from Equation[3] As 
the oscillation probability P(v^ — > v M ) depends on E v j cos 0, the natural choice for a variable in which the channel 
i can be defined is the ratio between a quantity which depends on the neutrino energy and the reconstructed zenith 
angle, ©«. As explained in Section |5J the energy-dependent variable is the observed muon range in the detector. The 
oscillation parameters are extracted by a x 1 minimisation which is detailed in Section|7] 



3. The ANTARES Detector 

A detailed description of the ANTARES detector can be found in The detector consists of 12 lines, equipped 
with photosensors, and a junction box which distributes the power and clock synchronization signals to the lines and 
collects the data. The junction box is connected to the shore by a 42 km electro-optical cable. The length of the 
detection lines is 450 m, of which the lowest 100 m are not instrumented. Their horizontal separation is about 65 m 
and they are arranged to form a regular octagon on the sea floor. They are connected to the junction box with the help 
of a submarine using wet-mateable connectors. Each line comprises 25 storeys each separated by a vertical distance 
of 14.5 m. The lines are kept taut by a buoy at the top of the line and an anchor on the seabed. The movement of the 
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line elements due to the sea currents is continuously measured by an acoustic calibration system with an accuracy of 
10 cm H. 

Each storey contains three 45° downward-looking 10" photomultiplier tubes (PMT) inside pressure resistant glass 
spheres - the optical modules @]. Some of the storeys contain supplementary calibration equipment such as acoustic 
hydrophones or optical beacons J3]. 

The signals of each photomultiplier are readout by two ASICs. The charges and arrival times of the PMT signals 
are digitised and stored for transfer to the shore station The time stamps are synchronised by a clock signal 
which is sent at regular intervals from the shore to all electronic cards. The overall time calibration is better than 
0.5 ns 10. Therefore the time resolution of the signal pulses is limited by the transit time spread of the photomultipliers 
(cr ~1.3 ns) ifioll and by chromatic dispersion for distant light sources. All data are sent to the shore station. With the 
observed optical background rate of 70 kHz per PMT at the single photon level this produces a data flow of several 
Gbit/s to the shore. In the shore station a PC farm performs a data filtering to reduce the data rate by at least a factor 
of 100 lllll . Several trigger algorithms are applied depending on the requested physics channel and on the observed 
optical noise. 



4. Data and Simulations 

The present analysis is based on data taken with the ANTARES detector between March 2007 and December 
2010. Until December 2007 ANTARES operated in a 5-line configuration, followed by several months of operation 
with 10 installed detector lines. The detector construction was completed in May 2008. All physics runs taken under 
normal conditions have been used. The events selected by two tight trigger conditions are used fl. The analysed 
sample consists of 293 million triggers, dominated by atmospheric muons, corresponding to a detector live time of 
863 days. 



Downgoing atmospheric muons were simulated with the program MUPAGE 112LI13I1 which provides parametrised 
muon bundles at the detector. Upgoing neutrinos were simulated according to the parametrisation of the atmospheric 
Vp flux from H 1 411 in the energy range from 10 GeV to 10 PeV. The Cherenkov light, produced inside or in the vicinity of 



the detector instrumented volume, was propagated taking into account light absorption and scattering in sea water 111511 . 
The angular acceptance, quantum efficiency and other characteristics of the PMTs were taken from |@1 and the overall 
geometry corresponded to the layout of the ANTARES detector J2l- The optical noise was simulated from counting 
rates observed in the data. At the same time, the definition of active and inactive channels has been applied from 
data runs as well. The generated statistics corresponds to an equivalent observation time of 100 years for atmospheric 
neutrinos and ten months for atmospheric muons. 



5. Reconstruction 

The algorithm used for the muon track reconstruction is described in flal . It assumes a simplified detector geom- 
etry composed of straight vertical detection lines. The method combines a strict selection of direct Cherenkov photon 
hits which are grouped around "hot spots" at each detector line with a^ 2 -like fitting procedure. A hot spot corresponds 
to a signal of at least 4 photoelectrons seen on two adjacent storeys of the same detector line within a narrow time 
window of less than 100 ns. Only hits on detector lines with such a hot spot are used in the track fitting. If the selected 
hits occur only on one detector line, a single-line fit is performed. No azimuth angle is determined in this case due to 
the rotational symmetry of the problem. This does not affect the present measurement, as the oscillation probability 
does not depend on the azimuth angle (see Equation 0. If selected hits occur instead on several detector lines, a 
multi-line fit is performed which provides both the zenith and azimuth angles of the track. The inclusion of single-line 
events is a special feature of the reconstruction method and allows to significantly lower the energy threshold of the 
final atmospheric neutrino sample. Whereas for multi-line events the threshold energy of the final neutrino sample is 
about 50 GeV due to the 65 m horizontal gap between lines, single-line events are reconstructed down to 20 GeV for 
nearly vertical tracks, accessing events with E v /L values close to the first oscillation maximum. 

The neutrino energy is estimated from the observed muon range in the detector. The selected hits are sorted 
according to their vertical position, z. The z-coordinates of the uppermost and lowermost hits, Zmax an d z m m, together 
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with the reconstructed zenith angle, Q«, allow to define an approximate muon range 



S = (Z,„ax ~ Zmin) I COS & R . 



(5) 



When considering the geometry of the Cherenkov light cone, Equation[5]is exact only for strictly vertical tracks. It has 
been verified that the use of a more sophisticated range definition does not improve the precision of the measurement 
for the oscillation parameters. S is used to estimate the muon energy by taking into account the ionisation energy loss 
of 0.2 GeV/m for minimum ionising muons in sea water 



E R = S • 0.2 GeV/m. 



(6) 



The derived energy from the visible muon range in the detector can be considered as a lower limit of the actual 
neutrino energy. The presence of a hadronic shower at the neutrino vertex is ignored, as well as the fact that the muon 
might leave or enter the detector, making only a fraction of its actual range available for measurement. The normalised 
difference between the true neutrino energy, E v , and Er is shown in Figure Q] for all selected simulated events (see 
Section |6]l with E v < 100 GeV. The mean value of 0.45 illustrates the fact that Er measures on average about half of 
the neutrino energy. The RMS of the distribution is 0.22. 
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Figure 1: Difference between the true neutrino energy, E v , and En normalised by E v for low energy events (E v < 100 GeV) from the final event 
sample of Section|6] 



6. Event Selection 

Downgoing atmospheric muons dominate the ANTARES event sample on the trigger level as illustrated in TableQ] 
Neutrinos contribute less than 10 here. In principle a simple cut in the reconstructed zenith angle should be enough 
to separate both classes. But as seen from TableQ] this is not sufficient. Atmospheric muons are often seen as bundles 
and they can be accompanied by hard stochastic processes such as bremsstrahlung. Both effects complicate their 
correct reconstruction and a certain fraction of downgoing atmospheric muon events are misreconstructed as upgoing. 
Some cuts on the quality of the reconstructed tracks are needed to reduce this contamination and derive reliably 
neutrino oscillation parameters from the data set. The goodness of the track fit is measured by the "normalised fit 
quality" as introduced in 11611 . a quantity equivalent to a^ 2 per number of degrees of freedom (NDF). The selection 
cuts, which are described below, have been obtained from a blind analysis. The single-line data sample has been kept 
blind, thereby masking a possible oscillation signature. 

For the multi-line selection, only events which have hits on more than 5 storeys are kept to allow a non-degenerate 
track fit. Further, the fit must not converge on a physical boundary of any of the fit parameters. As the contamination 
of misreconstructed atmospheric muons is particularly strong close to the horizon, a further condition, cos &r > 0.15, 
is imposed, i.e. tracks closer than 9° to the horizon are excluded. 
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Fit quality Fit quality 

Figure 2: Normalised fit quality of the final multi-line (left) and single-line (right) samples. Data with statistical errors (black) are compared to 
simulations from atmospheric neutrinos with oscillations assuming parameters from yj (red) and without oscillations (green) and atmospheric 
muons (blue). For a fit quality larger than 1.6 (multi-line) or 1.3 (single-line) the misreconstructed atmospheric muons dominate. The arrows 
indicate the chosen regions. 



The distribution of the normalised track fit quality of the resulting multi-line event sample for data and simulations 
is shown in Figure [2] (left). The neutrino Monte Carlo samples are scaled down by an overall normalisation factor 
r = 0.86 as obtained from the fit (see Section [8). This is well within the uncertainty of the atmospheric neutrino 
flux model flTtl . Figure [2] (left) shows that a cut on the normalised fit quality allows to cleanly separate the upgoing 
neutrinos from the downgoing muons. In order to have a contamination of misreconstructed atmospheric muons below 
5%, a cut value of 1.3 is chosen. 

For the single-line selection, events which have hits on more than 7 storeys are kept. This yields a minimal track 
length for a vertical upgoing muon of about 100 m, which can be produced by a muon of 20 GeV. Further cuts are 
identical to the multi-line selection. 

The distribution of the normalised track fit quality of the resulting single-line event sample for data and simulations 
is shown in Figure [2] (right). The neutrino Monte Carlo samples are again scaled down by a factor r = 0.86. Figure|2] 
(right) shows that also for this data set a cut on the normalised fit quality allows to cleanly separate the downgoing 
muons from the upgoing neutrinos. In order to have a contamination of misreconstructed atmospheric muons below 
5%, a cut value of 0.95 is chosen. 





Multi-line 


Single-line 




Data 


vMC 


//MC 


Data 


vMC 


^MC 


All 


1.42- 10 s 


8755 


1.23- 10 8 


1.51 ■ 10 8 


8242 


1.10- 10 8 


Nstorey > N c „, 


1.33 • 10 8 


8248 


1.18- 10 8 


4.44 • 10 7 


1260 


3.03 • 10 7 


Fit boundary 


1.32- 10 8 


8150 


1.17- 10 8 


4.31 • 10 7 


1242 


2.93 • 10 7 


cos0 s > 0.15 


2.74 • 10 6 


5512 


1.84- 10 6 


7.97 ■ 10 5 


1116 


6.96 ■ 10 5 


Fit quality cut 


1632+40 


1971+6 


52 ± 12 


494 ± 22 


651 ±3 


28 ±9 






1910 + 6 






557 ±3 





Table 1: Event reduction due to the cuts used. Statistical errors are given for the final data set. The effect of oscillations with parameters from Q 
is taken into account only for the values given in the very last row. 

The effect of the different selection cuts in the two channels is detailed in Table Q] for data and the Monte Carlo 
sets. Satisfactory agreement between data and Monte Carlo numbers is observed at all cut levels. 

Events from v e charged current (CC) interactions as well as neutral current (NC) interactions produce cascade-like 
event topologies. The spatial extension of these cascades, which are composed of the hadronic and electromagnetic 
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showers at the neutrino interaction vertex, does not exceed more than a few meters. This is in contrast to the lowest 
energetic selected muons, which have still a range of at least 100 m in sea water and therefore a distinct topology. 
The above described selection cuts, tuned to suppress the atmospheric muon background, reduce efficiently also the 
contribution from NC and v e interactions. Their contribution to the final event sample is estimated to be less than 
one event in the multi-line channel and 6 events in the single-line channel. Both numbers are significantly smaller 
than the residual background contribution from atmospheric muons. 

The zenith angle of the final neutrino sample is reconstructed with a precision of 0.8° for multi-line events and 
3.0° for single-line events (median of the angular error distribution with respect to the true neutrino direction). The 
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E Y /cos0 (GeV) 

Figure 3: Distribution of £,,/ cos & for the selected events of the atmospheric neutrino simulation. The solid lines are without neutrino oscillations, 
the dashed lines include oscillations assuming the best fit values reported in (3[ . The red histograms indicate the contribution of the single-line 
sample, in blue the multi-line events. 

distribution in E v j cos for the simulated selected atmospheric neutrino samples is shown in Figure[3] The expected 
event numbers are shown for the actual detector live time of 863 days. The effect of neutrino oscillations is seen by 
comparing the Monte Carlo curves with and without oscillations. The importance of the inclusion of the single-line 
events is evident. The expected deficit due to neutrino oscillations for E v j cos < 200 GeV amounts to 143 events, 
the majority of which (91 events) are in the single-line sample. These oscillated events are mainly converted into 
v T and might produce muons by a v T CC interaction followed by a decay r~ — > pTv^Vj (or charge conjugate). To 
estimate the appearance of muons from v T interactions in the final event sample, the energy dependent cross section 
ratio cr{v T CC) I criv^CC), about 0.5 for E Y = 25 GeV, the 17% branching ratio of the muonic r-decay and the soft 
energy spectrum of the resulting muons have to be considered. It is found that the expected deficit due to oscillations 
is reduced by less than 5 events. 

7. Systematic Uncertainties 

Systematic uncertainties lead to a correlated distortion of the Er/ cos®# distribution. They are implemented as 
pull factors in the x 1 function to be minimised: 

X 2 = J) [Ni ~ (1 + £)MC) L - (1 + 7,)MCf L ] 2 /erf + (e - t,) 2 /^. (7) 

i 

The sum extends over the bins in Er/ cos©« with N/ the observed number of events in bin i whereas MC\ L and 
MCf* L denote the expected event numbers from simulations in the single-line and multi-line channel respectively and 
cr, is the corresponding statistical error of bin i. The simulated event numbers depend on the oscillation parameters 
according to Equation|4] The quantities e and i] are two pull factors which allow to renormalise the two channels. The 
variations of e and 77 are constrained by cr R , which is introduced below. 



7 



A large class of uncertainties modify the overall normalisation, i.e. they act on e and rj in the same way. The 
normalisation of the neutrino flux model as well as the uncertainty on the absolute neutrino detection efficiency of the 
ANTARES detector fall into this class. Equation [7] does not constrain such uncertainties, as indicated by the absence 
of terms such as e 2 /<x 2 or ?7 2 /0" 2 , i.e. the relative overall normalisation between data and simulations is left totally free 
in this analysis. 

A second class of uncertainties might change the shape of the ErI cos©« distribution. Single-line events have 
typically lower E v j cos values (and thus correspondingly lower observed ErI cos &r values) than multi-line events 
(see Figure |3). A change in the ratio of the total number of events in each sample, Rq = N lL /N ML , would therefore 
reflect a shape change of the ErI cos 0s distribution. In terms of the pull factors this can be written as R(e,T]) = 
Ro(l + e)/(l + rj) as R()(l + e - 77). Therefore R(e, rj) varies with e-r\ if both pull factors are small. 

Several simulations have been performed with modified input parameters which may affect in a slightly different 



manner vertical and horizontal events (thereby acting on cos0r) as well as low and high energy events [18]. The 
average quantum efficiency of the phototubes was changed by +10% as well as their angular acceptance. The water 
absorption length of sea water was also varied by +10%. Further, the cuts in the normalised fit quality were varied in 
two steps of 0.05 around the chosen values thereby testing the stability of the analysis procedure. The spectral index 



of the atmospheric neutrino flux was varied by +0.03 as suggested in II 1911 . As a result a set of values of R was derived. 
It was observed that R remains stable within 5% of its original value, R$. This value <tr = 0.05 is used in Equation|7] 
to constrain relative variations of e and 77. This accounts for the maximal variation of the shape of the ErI cos ®r 
distribution due to the considered systematic effects. 

8. Results 

After applying the cuts defined in Section|6]the final event numbers obtained are given in TableQ] Figure 2] (left) 
shows the resulting ErI cos &r distribution for data and simulations. Figure |4] (right) shows the fraction of measured 
and simulated events with respect to the non-oscillation Monte Carlo hypothesis and indicates a clear event deficit for 
ErI cos 0« < 60 GeV, as expected assuming atmospheric neutrino oscillations. Restricting the y 1 minimisation to 




100 120 140 
E R /cos0 R (GeV) 




100 120 140 
E R /cos0 R (GeV) 



Figure 4: Left: Distribution of ErI cos Qr for selected events. Black crosses are data with statistical uncertainties, whereas the blue histogram 
shows simulations of atmospheric neutrinos without neutrino oscillations (scaled down by a factor 0.86) plus the residual background from atmo- 
spheric muons. The red histogram shows the result of the fit. Right: The fraction of events with respect to the non-oscillation hypothesis. Same 
color code as for the left figure. 

parameters in the physically allowed region yields the red curve of Figure |4]with Am 2 2 = 3.1 • 10~ 3 eV 2 and sin 2 2023 = 
1.00. The corresponding pull factors are e = -0.138 and 77 = -0.142. This is used as an overall normalisation 
parameter r = 0.86 in Figures|2]and|4] The fit is performed in 25 bins: the 24 bins shown in Figure|4]plus one overflow 
bin which contains 299 events with ErI cos©# > 140 GeV. The fit yields ,y 2 /NDF = 17.1/21. When imposing the 
world average oscillation parameters, Am 2 , = 2.43 ■ 10~ 3 eV 2 and sin 2 26*23 = 1 from yD, ^ 2 /NDF = 18.4/21 is 
found. 
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For the non-oscillation hypothesis, i.e. sin 2 2023 = 0, ^f 2 /NDF = 31.1 /23 is obtained. The pull factors in this case 
are e — -0.302 and rj = -0.196. The event deficit in the single-line channel is seen here as e becoming lower than rj. 
Requiring in addition e — r\ the x 1 increases further to ^ 2 /NDF = 40.0/24, which has a probability of only 2. 1 %. 

This measurement is converted into contours of the oscillation parameters and is shown in Figure |5] With 68% 



> 

CI 

C 6 

CM 
CM CO 

E 



1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



-6 



-5 



5- 



4 - 



3 - 



2- 



1 ' ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ "1 

t) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin 2 29 23 




Figure 5: 68% and 90% C.L. contours (solid and dashed red lines) of the neutrino oscillation parameters as derived from the fit of the ErI cos Or 
distribution. The best fit point is indicated by the triangle. The solid filled regions show results at 68% C.L. from K2K (2oll (green), MINOS l2lll 
(blue) and Super-Kamiokande f23l (magenta) for comparison. 

C.L. sin 2 2023 > 0.70 is found and Am 2 2 is constrained to values in the range [2.2,4.2] • 10~ 3 eV 2 . If maximal mixing 
is imposed (sin 2 2023 = 1), the obtained range of Am 2 , is 

Am 2 2 = (3.1 ±0.9)- 10~ 3 eV 2 . (8) 
The results are in agreement with other measurements from K2K (2oTl . MINOS f2lll and Super-Kamiokande f2"3l . 



9. Conclusions 

Based on data taken by the ANTARES neutrino telescope from 2007 to 2010, constraints on the neutrino os- 
cillation parameters sin 2 2023 and Am 2 2 have been derived. If maximal mixing is assumed, a value Am^ 2 = (3.1 + 
0.9) • 10~ 3 eV 2 is obtained. The result agrees well with current world data and demonstrates a good understanding of 
the performance of the ANTARES telescope at its lowest accessible energies. It is the first such measurement by a 
high energy neutrino telescope and underlines the potential of future low energy extensions of the existing neutrino 
telescopes for such physics. 
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